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Abstract 



We present a number of arguments relating magnetic disorder to center disorder, in pure 
Yang-Mills theory in D=3 and D=4 dimensions. In the case of the D=3 Georgi-Glashow 
model, we point out that the abelian field distribution is not adequatedly represented, at 
very large scales, by that of a monopole Coulomb gas. The onset of center disorder is 
associated with the breakdown of the Coulomb gas approximation; this scale is pushed off 
to infinity in the QED 3 limit of the 3D Georgi-Glashow model, but should approach the 
color-screening length in the pure Yang-Mills limit. 
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1 Introduction 



The Zn center of an SU(N) gauge group is associated with the confinement properties of 
of a pure gauge theory in a number of ways. It is well known that the finite temperature 
confinement /deconfinement transition can be regarded as the breaking of a global 
symmetry in a volume with a compactified time direction. In addition, as shown by 
't Hooft the VEV of a Zn vortex creation operator can be interpreted as an order 
parameter for confinement, dual to Wilson loops, in SU(N) gauge theory. It was also 
suggested many years ago that "thick" Zn vortices are responsible for the area-law falloff 
of Wilson loops |l|-f§, and recently there have been a number of numerical investigations 
which support this idea [|IO|-|T2"lj . 

The notion that confining ("magnetic") disorder is center disorder may also be sup- 
ported by some simple observations, presented in section 2, regarding the behavior of 
holonomy probability distributions in Yang-Mills theory. We point out that the holon- 
omy distribution approaches a random distribution on the group manifold as loop size 
increases; however, the approach to the random distribution is far more rapid among the 
center elements than among elements of the coset. We also show, with the help of the 
lattice strong-coupling expansion, that while center elements within a large area fluctuate 
independently, this is not true of fluctuations in the coset for D > 2 dimensions. 

The 3D Georgi-Glashow model (GG3) is interesting in this context for several reasons. 
The confinement mechanism in this theory is believed to be essentially that of compact 
QED 3 , at least in some region of the coupling parameters; one therefore expects that 
confining disorder is U(l) disorder. There is only one phase in GG3, but there are two 
special limits: Compact QED3 is obtained in the limit where the mass of the W-boson 
becomes infinite, while pure Yang-Mills theory (YM3) is obtained in the limit where the 
adjoint scalar effectively decouples from the gauge field. Since the 3D Georgi-Glashow 
model interpolates smoothly between QED 3 and YM 3 , a natural question to ask is what 
happens to center disorder in GG3 as we move away in parameter space from the pure 
Yang-Mills limit. 

This question is taken up in section 3, where we point out a qualitative difference 
between confinement in compact QED3, and confinement in the 3D Georgi-Glashow model 
at large scales. In the case of compact QED%, we show via saddlepoint methods that 
double-charged loops have twice the string-tension of single-charged loops, while for GG3, 
the double (abelian) charged loops must ultimately be screened by massive W-bosons. 
As a consequence, the effective abelian theory corresponding to the 3D Georgi-Glashow 
model, obtained after integrating out the charged bosons, is not adequately represented 
by a Coulomb gas of 't Hooft-Polyakov monopoles. Our main point is that the massive 
W-bosons of GG3 are not just spectators whose effect on vacuum fluctuations, beyond the 
range , is negligible; in fact the W-bosons must strongly affect the vacuum distribution 
of abelian flux at large distance scales. At these large scales, it appears that confining 
disorder in GG3, as in the pure Yang-Mills theory, is associated with Z 2 (rather than U(l)) 
disorder. In the QED3 limit, the onset of Z2 disorder is pushed off to infinity, while in the 
pure Yang- Mills limit, it roughly coincides with the onset color screening. 
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2 Confining Disorder as Center Disorder 



A Wilson loop is understood as measuring the response of the vacuum to the introduction of 
heavy sources, but it can also be viewed as providing information about field fluctuations in 
the ground state, in the absence of external charges. Consider, in particular, a gauge theory 
with matter fields in the fundamental representation. The asymptotic perimeter-law falloff 
of the Wilson loop is explained by the binding of matter quanta to the external charge, 
forming a color singlet. On the other hand, imagine integrating out the matter fields, 
leaving an effective action involving only the gauge fields. It is then clear that the effect of 
the virtual matter fields is to modify the probability distribution of gauge-field fluctuations, 
such that confining configurations, which would normally induce an asymptotic area-law 
falloff in the loop, are suppressed. 

In discussing the probability distribution of gauge fields, with or without the presence of 
matter fields, it will be helpful to introduce a gauge-invariant operator which is somewhat 
more general than a Wilson loop. Consider, for simplicity, a lattice pure-gauge theory with 
an SU(2) gauge group, and let U(C) denote the path-ordered product of link variables along 
a closed loop C (the holonomy). The expectation value 

P c [g] = <5[g,U]>=<J2xMxAU(C)]> 

j 

= l + E^(C)XiM (!) 

is the probability density, on the SU(2) group manifold, that the loop U{C) equals the group 
element g. The sum over j runs over group representations, the Xj(fiO are SU(2) group 
characters, and Wj(C) =< Xj[U(C)] > is the VEV of the Wilson loop in representation j. 
As loop C becomes large, this holonomy probability distribution approaches the random 
distribution, Pc(g) — ► 1, and it will be useful to focus on the deviation of Pc, denoted 
Pc[g], from the random distribution 

Pc[g] = Pc[g] - i (2) 

Since Pc[g] is gauge-invariant, it can only depend on the eigenvalues of the unitary matrix 
g, and Pc has flat directions on the group manifold corresponding to g — > g' — ugv) '. In 
the particular case of SU(2), Pc only depends on Tr(g), and we will be interested in how 
this dependence fades away as the loop C becomes large. 

A Wilson loop Wj(C) can be thought of as a moment of the probability distribution 
Pc(g)- It is expected that, in D = 3 and D = 4 dimensions, planar Wilson loops have the 
asymptotic form 

W (C) = I e M~ a AC) - PjPiC) - Cj\ J = half-integer 
^ ' ) exp[—fj,jV(C) — Cj] j = integer ^ ' 

at any lattice coupling. Note that all half-integer representations have the same string 
tension, and all integer representations have zero string tension. This is a well-known con- 
sequence of color-screening, which seems (from numerical studies) to set in somewhat after 
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the onset of confining behavior. The perimeter term reflects both short-range, perturbative 
contributions, roughly proportional to the quadratic Casimir j(J + 1), and also, for j > |, 
the bound-state energy of gluons required to screen the color charge to its minimum value 
(either j — or j — The constant Cj, which increases with j, can be attributed to the 
rapid initial falloff of higher-representation loops in the so-called "Casimir-scaling" region, 
before the onset of color-screening [[H|,[14] . At the point where the asymptotic behavior sets 
in, the higher-representation loops have already fallen to a rather small value as compared 
to lower- representation loops, and this fact is accounted for in the constant Cj. Since both 
fij and Cj increase with j, it follows that for large loops 

W 1/2 (C) » W 3/2 (C) » W 5/2 (C) » ... (4) 

and 

W X {C) > W 2 {C) > W 3 (C) > ... (5) 
It also follows, for sufficiently large loops, that 

Wi(C)>W 1/2 (C) (6) 

since the rhs falls off asymptotically with area-law behavior, and the lhs only falls off with 
the perimeter law.Q 

Using we have the leading behavior 

Pc(g)*Xi(9)ew[-»iV(C)- Cl } (7) 

and the approach to the purely random distribution follows a perimeter-law, rather than 
the area law which might have been expected. In contrast, in D=2 dimensions 

Wf{C) = (2j + 1) expi-ajAiC)] 

o 3 = - lo s^y- (8) 

and therefore, asymptotically, 

P^ D (g)^Xi/2(g)ex V [-a 1 _A(C)] (9) 

This 2D distribution, unlike the 3D and 4D distributions, approachs the random value via 
an area-law falloff. 

There is, however, a hidden area-law approach to randomness also in the D = 3, 4 
group distributions. Let us extract a center element from the holonomies 

Z[U{C)\ = signTr[[/(C)] (10) 



1 It should be noted that condition (||) has yet to be verified numerically, at least at zero temperature. 
The onset of color screening appears to be at the edge, or perhaps just beyond, the range of current 
numerical simulations (cf. Michael in ref. [Q). The condition can be verified at small 0, using the lattice 
strong-coupling expansion. 
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and ask for the probability Pc{ z ) that z = Z\U{C)\, where z = ±1. This is given by 



Pc(z) =< ~(1 + z x signTr[[/(C)]) > (11) 

The pure-random value is Pc{ z ) — 1/2, and again we remove this constant to define the 
deviation Pc(±) = -Pc(i) — § from pure-random. Then, from the character expansion 

signTr[f/(C)] = £ a iXi [[/(C)] 

dj = J dg signTi[g]x j(g) (12) 

and again applying (f|), we find 

P c {z) « z|- exp[-a 1/2 „4(C) - // 1/2 P(C) - c 1/2 ] (13) 

The conclusion is that, although the overall holonomy probability distribution Pc[g] 
approachs the random value via a perimeter falloff in D=3 and D=4 dimensions, the 
probability that Tr[C/(C] has one or the other sign approachs the random distribution 
via an area-law falloff. In D=2 dimensions there is no such distinction between Pc{g) 
and Pc{z); both probabilities have an area-law falloff. The strong implication is that 
fluctuations in the center element, which distinguishes between two cosets of the group 
characterized by the sign of Tr(g), are the fluctuations characteristic of confining disorder 
in D=3 and D=4 dimensions. To go further, however, we will need to resort to the lattice 
strong-coupling expansion. 

We note in passing that one finds the relation 

< signTr[f/(C)] >^ ^W 1/2 (C) (14) 



from the character expansion (|1"2"D, as a consequence of the inequalities (f4|). This explains 
the rather mysterious equality of potentials extracted (i) from the Wilson loops Wi/2(C); 
and (ii) from the sign of Wilson loops < signTr[[/(C)] >, which was found recently in 
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2.1 Z 2 Disorder at Strong Coupling 

It is often said that confinement in strong-coupling lattice gauge theory is simply a matter 
of plaquette disorder: Group elements associated with loops around nearby areas (the 
plaquettes) fluctuate independently, leading to an area law falloff for the Wilson loops. 
This is certainly a correct statement of the situation in D = 2 dimensions. However, as we 
will now show, there are some important qualifications to be made in higher dimensions. 

Let us consider a very large planar loop C whose minimal area A(C) is subdivided into 
some number n of subareas A(Ci), encircled by loops Ci (Fig. [I]). If all loops Ci are large, 
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Figure 1: Subdivison of a large loop C into smaller regions bounded by loops {Cj}. Note 
that the sum of interior perimeters may be much greater than the total perimeter of C. 



and the coupling is strong, the question is to what extent the holonomies U{Cj) fluctuate 
independently. Let F[g] denote any class function with character expansion 

/••..'/ X/A,(.</) (15) 

The test for whether the U (Cj) fluctuate independently is whether or not the VEV of the 
product of F[U(Ci)) equals the product of the VEVs, i.e. 

<nF[U(C t )]>=Yl<FiU(C t )}> (16) 

i i 

In D=2 dimensions, it is easy to verify that the equality holds 

<f[F[U(C i )]> « (n2/i/ 2 )exp[-a 1/2 „4(C/ 

n 

= H<F[U(C l )]> (17) 



1=1 



i=l 



so the group elements U{Ci) do seem to fluctuate independently in each subregion. 

In D > 2 dimensions the answer is different. Suppose that each Cj satisfies 
A(Ci) ^> 2V(Ci) in lattice units. In that case, we find from the strong-coupling expansion 
the leading contribution to the VEV of products^] 

n / 1 \ n— 1 

<IIWi)]> « U) K<xx[u{c)]> 

ftexp[-VLV{C)] (18) 



1 \ n ~ ^ 



2 We neglect here certain sub-leading, shape-dependent terms in the exponent 
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where /i = 4cri/2. The powers of 1/3 are due to the fact that this contribution is highly 
non-planar, and would vanish in the large-N limit. The planar contribution, however, has 
an area law falloff, and for any finite N it is negligible for large loops. Likewise, for the 
product of VEVs, we have the leading contribution 



n<^(c*)]> « n/i<xi[^(c i )]> 



i=l 



= /fexp 

However, for n ^> 1, 



(19) 



(20) 

i=i 

from which we see that the exponential falloff of < Jlj F[U(Ci)] > and Ili < F[U(Ci)] > is 
quite different, already at the leading terms in the exponents. The conclusion is that the 
group elements U(Ci) do not, in fact, fluctuate independently; the deviations from pure 
random in the set of sub-loop probability distributions must be correlated. 

On the other hand, since the area-law falloff of Wilson loops is supposed to be due 
to magnetic disorder; some component of magnetic flux should be fluctuating (nearly) 
independently. It is easy to check, at strong coupling, that the center elements Z[U(Ci)] 
have the required property. The leading behavior for the VEV of the product is given by 

<l[Z[U(C i )> = <nsignTr[C/(Ci)] > 

i i 

- <n(^)xi/ 2 [^)]> 

/ 3 \ 11 

« e W [-a 1/2 A(C)} (— j (21) 

Similarly, for the product of the VEVs 

H<Z[U(C i )> = n < signTr[?7(a)] > 

i i 

« IIl^J^ph^iMCa)] (22) 



i=l 



Comparing ( |2TD and (0), we see that 

<J\Z[U{C i )]>nJ[<Z[U{C i )]> (23) 

i i 

The center elements of the holonomies therefore fluctuate independently. 

Any SU(2) class function F[U(Ci)] can be expressed as a function of of the sign and 
the modulus of Tr[C/(Cj)] as follows: 

F[U{Ci)\ = F^iUiQ)}] + Z[U(C t )]F 2 [Ti 2 {U(C t )}} (24) 



7 



VEVs of products of F\ and F 2 do not factorize, as seen in eqs. (pT8|-pO|), while products of 
the Z\U{Cj)\ do factorize, as seen above in (|23|) . The conclusion is that there is magnetic 
disorder in the center elements Z of the SU(2) group holonomies, but not in the coset 
elements, which depend only on Tr 2 [[/(C)]. 

All of this has some bearing on the question of what are the relevant confining con- 
figurations in non-abelian gauge theory. Confining configurations, whatever they may be 
in D=3 and D=4 dimensions, must have the property of disordering the signs of Wilson 
loops, but not disordering the absolute values of those loops. At least, we have seen that 
this must be true at strong coupling. A class of configurations with these properties is the 
"spaghetti vacuum" of center vortices, proposed twenty years ago, in various forms, by 't 
Hooft [^|], Mack 0, the Copenhagen group @-f|, and others 0-0. A center vortex, linking 
loop C, has the property of sending U(C) — > zU(C); such configurations can only disorder 
the center elements (i.e. the signs of Wilson loops, in SU(2)), leaving the rest of the group 
distribution untouched. This seems to be exactly what is needed. 

On the other hand, it is not at all excluded that there could be some other degrees of 
freedom associated with loops that also fluctuate independently, and contribute to magnetic 
disorder. In our discussion so far we have only considered the holonomies U(C), but we 
could also imagine using, e.g., an adjoint Higgs field (either elementary or composite) to 
construct other types of loop elements that might fluctuate independently. In fact, this 
is the general idea behind monopole confinement: An adjoint Higgs field is used to single 
out a U(l) subgroup of SU(2), and it is U(l) group elements, associated with loops Cj, 
that are disordered via a monopole condensate. From this point of view, the Z 2 disorder 
could be just a subset of a more general U(l) disorder. The simplest and most explicit 
proposal for monopole confinement (and U(l) disorder) in a non-abelian gauge theory is 
due to Polyakov, in his analysis of the D=3 Georgi-Glashow model JT5 , IB]. Since GG3 



interpolates smoothly between QED 3 and pure YM%, we would now like consider if there 
is any remnant of Z 2 disorder in the 3D Georgi-Glashow model. 

3 Double-Charged Loops in GG3 

The full lattice action of the Georgi-Glashow model is a function of the gauge field variables 
Up(x) and the adjoint Higgs field variables <p{x) 

s[u,<p] = ±PgY,t±[uuuW] 

plaq 

+ ±p H ^Tt[Up(x)<l>(x)Ut(x)<p(x + [i)] 

~ eI^^I+^^^I- 1 ) 2 } ( 25 ) 

where 

</>(x) = y £r(x)<T a (26) 



a=l 



s 



The naive continuum limit is obtained from the lattice action by the following scaling: 



U,(x) = e iaA ^\ (A^^A^), <p = a 1 / 2 (3 1 R /2 ¥ >, (27) 
and the identification of continuum couplings g, A and fi by 

= — , — = -Aa, = -- /ia . 28 

g 2 a p H 4 f3 H 2 



Thus the continuum action becomes^ 



1, 



/ 2 



^i+^(M 2 4^ 2+ i A ^ 4 ]' (29) 



where -D M yj = + i[A^, <p]. It is obtained for 

Pg^oo, (3r^0, Ph^^ (30) 

the approach monitored by the lattice spacing a. A more precise description is obtained 
by taking into account that although the coupling constants are not renormalized in three 
dimensions, the masses M\y and M v are renormalized. We refer to |18| for formulas 
which include the one-loop mass renormalizations. The tree-formulas describe correctly the 
qualitative features which have our interest. In the "broken" region, where the approach 
of (3h to 1/3 is from above, it corresponds mass M w for the charged vector particle and a 
mass M v for the neutral scalar, given by: 

M 2 w = g 2 V l M 2 = 2^ 2 C , & = £ (31) 

If we in the limit fl30l) increase (3h away from 1/3, it corresponds formally to a "contin- 
uum" limit where 

M w ~ - M v ~ -, (32) 
a a 

i.e. a limit where only the free photon field is present as a physical excitation. More 
generally, taking j3 G and j3 R fixed and (3h — > oo have the same effect, except that the 
resulting effective theory will be compact lattice U(l), where we obtain the continuum free 
photon field only in the limit (3g oo. For future reference let us note that the tree- value 
formulas (p8|) and ( pi]) lead to the following expressions for dimensionless ratio 



9 2 4 I 20 R 



(33) 



The actual region above (3h — 1/3 where Mw is not of the order of the inverse lattice 
spacing a is very narrow, in agreement with numerical simulations [ITS . 



J A proposal for the effective action in dual variables is found in ref. p7| . 
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Going to the unitary gauge, we write 



S ug [U, p] = S[U, </>(x) = p(x)a 3 ] (34) 

The unitary gauge action S ug still has a residual £7(1) gauge symmetry. We may factor the 
SU(2) link variable into a matrix which transforms under the residual symmetry like a 
£7(1) gauge field, and a piece C M transforming like a double-charged matter field [19| 



cos ipe l9 sin ipe lx 
- sin ipe~ lx cos (pe~ ld 



cosy? sin^e* 7 \ ( e ld 
— sin (pe~ 11 cos I 1 e - *' 



with lattice measure 



VI - cc* c \ e ie 
-c* VI - cc* J 1 e" ie 



(if/ = — - / dip cos 99 sin tp / ^6* / (36) 
2ti 2 Jo J -tv J -tt 

The effective abelian action S^/LA] for GG 3 is then defined by 

e Se ff [A] = I Y[dp(x) Y[dip^x)d Xfl {x) cos^(x) sin^(x) e s -^ (37) 

The question we wish to raise is whether the Euclidean quantum theory of S" e //[y4], at 
large scales, is correctly represented by a monopole Coulomb gas, as proposed in [[L5[], since 
this is the generic model of monopole confinement. In particular, consider the "abelian" 
Wilson loops 

A n {C) = Tr[{AAA...A) n ] (38) 

corresponding to closed loops of heavy particles carrying n units of the abelian electric 
charge, with 

< A n (C) >= i J DA A n (C)e s <ffW (39) 

The crucial point is that in GG3, the string-tension a n of loops < A n {C) > should vanish 
asymptotically, if n is an even integer. The reason is simply that the Georgi-Glashow model 
contains massive W-bosons which carry two units of electric charge; these correspond to 
the Cfj,(x) degrees of freedom. The W-bosons are able to screen static sources carrying an 
even number of unit electric charges. If a flux tube were to form between n=even charged 
sources, with some string tension a n = T then at a separation of roughly 

nMw , x 

L = -f- (40) 

charge- screening by W-bosons becomes energetically favorable and the flux-tube breaks, 
so a n = asymptotically. Although S^z/LA] contains only the abelian gauge field, it must 



C^x)A^x) (35) 
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somehow incorporate this non-confinement of even charges, since we can always rewrite eq. 



as 



< A n (C) >= J DAD<pD x Dp A n (C)e s ^ [c ^ p] (41) 
where the charged fields are included explicitly. 

3.1 Strong coupling expansion 

While charge- screening may be expected in GG 3 on very general grounds, it is also possible 
to verify the effect explicitly in a strong-coupling expansion. Take, for simplicity, (3r = 00 
so that the modulus of the Higgs field is "frozen" to p = 1, with [3h chosen large enough such 
that M\y 3> 1 in lattice units, and (3g small enough to allow a strong-coupling expansion. 
Expanding the action to 2nd order in c, c*, one easily finds for large, double-charged loops 
the perimeter falloff expression 

< A 2 {C) >= exp[-nV{C)] (GG 3 ) (42) 
with perimeter coefficient, extracted from the leading diagram, 

% 



H = - log 



/% 1 - 2(3 H e~ 2 ^ - e- 2 ^ 
8(3 H 1 - e- 2 ^ 



(43) 



In strong-coupling compact QED, of course, the answer is different. There are no charged 
bosons to screen the double-charged loop, and its value, for the Wilson action, is 

< A 2 (C) >= exp[-2aA(C)} (QED 3 ) (44) 

with a the string tension of the single- charged loop. For multiply charged loops < A n (C) > 
in strongly-coupled compact QED 3 , the string tension is in general n times the string 
tension for single- charged loops. Thus we have a qualitative difference, at least in strong- 
coupling, between compact QED 3 and S e ff[A] of the 3D Georgi-Glashow model, because 
in the latter effective abelian theory, there is no string tension for < A n (C) > when n is 
even. 

For compact QED 3 , it is possible to derive the result ( pBD also at weak couplings, 
which we do in the next section. This is an interesting result in its own right, since the 
existence of the n = 2 string tension in compact QED 3 has been questioned (ref. ||16|| , p. 
80), while, in numerical simulations, a finite value equal to twice the n = 1 value has been 
measured |2(J. 



3.2 Double-Charged Loops in a Monopole Gas 



It is well known that compact U(l) on a lattice can be written in a monopole gas repre- 
sention. If we use the Villain version of compact U(l) one obtains 



E 



exp 



m(r)= 



2vr2 
g 2 a 



/]m(r')G(r — r')m(r) 



(45) 
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where m(r) is an integer valued monopole field at the (dual) lattice size r, G{r — r') is the 
lattice Coulomb propagator in three dimensions, i.e. A 2 G(r — r') = —8 rr r, and g 2 a, a being 
the lattice spacing, is the temperature 1/(3 in the usual thermodynamic interpretation. 
One can represent the propagator G(r — r') by a Gaussian functional integral: 



n^xO) exp[-|-^^-(A M xW) 1 



x ex P 

m(r)=— oo 



2tt2 
g 2 a 



G(0) V] m 2 (r) + z ^ m(r)x(? 



(46) 



In the weak coupling limit (low temperature limit) where g 2 a — > we need only to maintain 
the first terms |m(r)| < 1 in the sum and we get 



7 ~ 

t-'mon ~ 



J]^dx(r) exp 



£ -(A M x) 2 + M 2 (l-cos X (r)) 



(47) 



4tt 2 V V 2 

where Mq comes from the propagator G(r — r') for coinciding arguments, (G(0) ~ 0.253 

2 



H) 



M 2 



87T 

— exp 
g*a 



2tt^ 



G(0) 



(4£ 



and if we interpret the monopoles as instantons, 2ir 2 G(0)/g 2 a can be viewed as the action 
of the instanton. 

Let us for notational simplicity carry out the following discussion in a continuum no- 
tation. All manipulations done in the following have a precise lattice analogy (see |21|). 
The translation is 

A^ad,, a 3 ]T^ / d\ (49) 

and we end up with 



7 ~ 

t-'mon ~ 



T>x(r) exp 



4vr 2 



^ (-(^ X ) 2 + M 2 (l-cos X (r))) 



where 



M 2 = a~ 2 Ml 



(50) 



(51) 



In the Coulomb gas approximation we are effectively integrating over the electromag- 
netic fields carried by the monopoles. They can be found from the monopole density m(r) 
by integrating 



d tl H t ,{r) = 2irm(r), i.e. H^r) = I J d 3 r' r~ ^rf ™.{r') 



(52) 



In particular we find, if C denotes a closed curve and S(C) a surface with C as boundary, 
that 



c 



dr^A^r) 



S(C) 



dS^r) H,(r) 



d 3 r Vs(C)(r) m(r), 



(53) 
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where 

^>M = -^/ S(C / s *M^r < 54) 

It is seen that rj(r) has the interpretation as a dipole sheet on the surface S(C). Thus, if 
C is a planar curve in the (x, y) plane and S(C) the planar surface with C as boundary 
curve, we have 

-d 2 r ]s{c) = 27T5'(z)6 s (x,y), (55) 

where 9s(c)( x , y) is one f° r a point inside the boundary C and zero for a point outside the 
boundary C. Close to the surface S we have 

r)(r) = 7rsign2 9 s (x,y), (56) 

i.e. it jumps by 2tt when passing the dipole sheet. 

We can now calculate the expectation value of a planar Wilson loop which carries n 
units of charge. From (|53| ) we obtain 

(A n (C)) = (e in $ dr ' iA ''W) = (^ e in J d3 ^S(C)(r)m(r)^ ^ 



where the expectation value is calculated with respect to the partition function (j45|) . Per- 
forming the same transformations which lead from ( f4"5|) to fl50|), the last term in ( |5TD leads 
to a translation of xi r ) such that one obtains 



,2 



(MQ) = J—j V x(r) exp[-f- / d 3 r (~{d»{ X - n Vs{c) ) 2 + M 2 (l - cos X (r)))] , (58) 

Since we consider the weak coupling regime the dominant contribution to the expecta- 
tion value (|58|) comes from the classical solution to the effective action in fl5"8"). In case we 
choose the planar loop to lie in the (x, y) plane we obtain from (|55|) : 

d 2 x = 2im5'{z)6 s (x,y) + M 2 sm X . (59) 

A solution to the homogeneous equation is (far from the boundary, suppressing a trivial 
(x, y) dependence): 

X (0 \z) =4arctane- M *. (60) 

Note that for z < it can be written as — 4arctane M2: + 2ir. Thus, for n — 1, eq. ( |59"D has, 
far away from the boundary C, the solution 

Xda = sign z ■ 4arctan(e~ M|z| )# s (:r,?/). (61) 

The important property of the solution (31) is that Xdii z ) ~* f° r \ z \ ~* 00 • This implies 
that it can be joined to the trivial solution \ = for |r| — > oo in i? 3 . 

Since the solution Xdi is given in terms of elementary functions one can calculate the 
corresponding action in (jj8|): 

d^x(^-(d^(Xda~ Vs(C)) 2 + M 2 (l — cosXdli 00 ))) = 8M Area(S) + perimeter contributions. 
2 (62) 
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In this way one obtains the famous area law of Wilson loops in the three dimensional 
Coulomb gas of monopoles since it follows from d58|) that 

r Q 2 i 
(A(C)) ^exp[-^8Marea(S)J. (63) 

Let us now consider the situation for double charged Wilson loops. In order to find the 
minimum of the action we should solve (]S9|) for n — 2. One could be misled to suggest the 
simple solution (far from the boundary of C) 

Xdl( z ) = 27rsign2;, i.e. x( z ) = 2 v( z )- (64) 

This clearly gives energy zero in the interior of the Wilson loop and seemingly no area law 
(as for the full Georgi-Glashow model). However, since this solution does not go to zero 
far from the Wilson loop, contrary to the solution ( |61]) for n = 1, we have to interpolate 

(2) (2) 

between the limiting values X c ia( z ) = ^ n an d Xda( z ) = m some region in space far 
away from the Wilson loop. This will cost an energy which is easily seen to be proportional 
to the area of the "domain wall" where the interpolation takes place. Thus the optimal 
situation is also here one where Xda( z ) ~* f° r \ z \ ~^ 00 • With this requirement it follows 
that we have to solve eq. fl59| ) with the boundary conditions that Xdl( z ) ~ * f° r N — >• oo 
and Xdl( z ) ~ ^v( z ) is differentiable for z = when passing the sheet. 

There is no such solution. But we can find approximate solutions with energies above, 
but arbitrary close to twice the energy corresponding Xdi- From (]5PD, (pT|) and ( |B"UD it 
follows that for z 3> 1/M 

X {2 \z) = 6{z) x i0 \z - z ) + d{-z){ X {0 \z + z ) - 2n) (65) 

is a solution to fl59|) with n = 2 except for exponentially small corrections, and its energy 
is twice that of Xdi except for exponentially small corrections. Further we see that this 
approximate solution behaves like (^) for \z\ « z . 

It is interesting that one can find a different kind of solution with the same features, 
namely that the energy of the classical solution can be arbitrary close to twice that of Xdii 
but never reach it. We have a free choice for the surface S(C), except for the requirement 
that C is the boundary of S. In particular, we could for the doubled charged Wilson loop 
choose two sheets separated a distance 2d and located in the z = ±<i-planes, except close 
to the boundary C. For each sheet we now have a discontinuity corresponding to rj(z). If 
d » 1/M it is clear that the solution, except for exponentially small corrrections, has to 
be 

x W(z) = e(-z) X ( 2( z + d) + e(z) X ( dl(z - d)W. (66) 

The energy becomes minimal (and equal two times that corresponding to Xdli z )) m the 
limit d — > oo. 

We conclude that the string tension for a double-charged Wilson loop will be twice the 
string tension of a single charged Wilson loop if we restrict ourselves to the Coulomb gas 
approximaton of functional integral. Clearly the arguments can be extended to n-charged 
Wilson loops. 
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3.3 Weak coupling limit of lattice GG3 

The monopole gas calculation above was a weak coupling expansion in the sense that g 2 a 
had to be considered small in order to make the truncation ( [47| ) . In particular this implied 
that the monopole "action" 2ir 2 G(0)/ g 2 a is large and the density of monopoles, 

2vr 2 G(0) 

p ~ exp' 



(-^), (67) 



9 a 

exponentially small. In the naive continuum limit, as defined by (f28"|)-(|31|), we can make 
contact to the similar instanton calculations in the continuum GG3 model. In that case 
the instanton (monopole) action is given by 

Smon = ^fe(\/g 2 ), (68) 
9 

where e(x) is a slowly varying function of x (e(0) = 47r). (It is seen that one obtains the 
QED 3 formulas in the limit where M\y ~ 1/a, as expected). A dilute instanton calculation 
is valid if the density of instantons, 

p ~ exp(-S mcm ), (69) 

is exponentially small relatively to the extension of the instantons (which is ~ 1/M^y). 
Thus the calculation in the last subsection is valid in GG 3 provided M w /g 2 ^> 1. One 
obtains a string tension for an n-charged Wilson loop 

T n ~ n e~ Smon x [subleading corrections] . (70) 

To the extent one can trust the tree-value formulas, the lattice inequality corresponding 
to Mw/g 2 3> 1 can be obtained from 



M w p H p G f3-(l-2p R )/p H 



2 



9 2 4 V 2/3 



R 



> 1. (71) 



This formula is most reliable for flc large (and [3r is small) and (3h close to (1 — 2/3 R )/3. 

We have seen above (see eq. ([40])) that we expect a perimeter law in GG3 for n-charged 
Wilson loops, n even, provided the linear extension L of the loop satisfies 

L > — ~ e s mo n x [ su bleading corrections] (72) 

J- n 

This length is much larger than the length scale 

ia = o~~ 1/2 ~ exp(^S mon y (73) 



set by the string tension, and it is also much greater than the average distance 

3 



R=p- 1 / 3 ~exv(ls mon ), (74) 
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between the monopoles. Typically we will have to go to distances larger than £ CT if we want 
to measure the string tension. However, for Mw/g 2 > 1 we have to move out exponentially 
many units of length £ ff before an n-charged string, n even, breaks: 

L 1 

— ~ exp(-S mon ) x [subleading corrections]. (75) 

In the lattice GG 3 model (|25|) the parameter (3 H allows us to interpolate continuously 
between compact QED 3 (large and pure Yangs-Mills theory {(3 H — > 0). From the 
tree-formula (|7T|) we see (for large /3o) how large (5h corresponds to a value M w /g 2 ^> 1, 
while Mw/g 2 drops to zero (in the tree-approximation) for /3# = (1 — 2/3r)/3. Below this 
value of (3h we have the "unbroken" coupling region of the Yang-Mills-Higgs system, where 
we expect the center Z 2 to play the dominant role in confinement and the monopole gas 
description is not valid at all. 



4 Z 2 Disorder in GG 3 

We now return to the question of Z 2 disorder in GG 3 . Considering only the abelian 
magnetic flux probed by loops A n (C), we can ask if disorder is distributed evenly in the 
U(l) group, or if it is only present in some subset of the degrees of freedom. Our procedure 
is the same as in section 3. Defining again the holonomy distributions on the compact U(l) 
group 

00 

P c (e iu >) = < 5[u, 9(C)] >=< ]T e in[ "- e(c)] > 

n=— 00 

Pc(0 = Pc(0-1 (76) 

we see that for compact QED 3 the approach to a pure random distribution has an area-law 
falloff 

P c (e^) ~ cos(^)e- CT - 4(c) (QED 3 ) (77) 
while in GG 3 , the approach goes, asymptotically, as a perimeter-law falloff 

P c (e^) ~ cos(2cj)e-^ {c,) (GG 3 ) (78) 

due to the different behavior of the n =even charged loops. However, once again, there is 
a hidden area-law approach to randomness also in GG 3 , since if we define 

Z\A[C)\ = signcos^C)] (79) 

and the probabilities 

P c {z) =< -(1 + zx sign[cos#(C)]) > (80) 
with z = ±1 and P c {z) = P c {z) - \, then in GG 3 

Pciz) ~ ze- aA{c) (81) 
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Turning to lattice strong coupling we again find, in complete analogy to the pure-gauge 
theory in section 3, that at leading order for compact QED 3 

< n F[A(Ci)] >^ II < F[A{Ci)\ > (QED 3 ) (82) 

i i 

while in GG 3 

<UF[A(C t ))> ~ 

i 

IJ < F[MCi)\ > ~ e -A*E 4 m) ( 83 ) 

It follows that, in contrast to QED 3 , the abelian loop elements do not fluctuate indepen- 
dently in GG 3 , even at very strong lattice coupling. On the other hand, the Z 2 elements 

<l[Z[A(C i )}> = IJe-S*^ 

i i 

= J[<Z[A{C i )\> (84) 

i 

do fluctuate independently in GG 3 , at strong coupling. 

The conclusion is that even in GG 3 , long-range disorder seems to be associated with 
a Z 2 , rather than a U(l), subgroup; there is again disorder in the sign, but not in the 
modulus, of loop elements cos 8(C). The inclusion of an adjoint Higgs field does not seem 
to change the fact that disorder, at large scales, is essentially a property of the gauge group 
center. In the last section we gave a qualitative description of the length scales in GG 3 
beyond which the Coulomb gas picture breaks down and where (as we have now argued) 
the magnetic disorder is center disorder. 



4.1 Extension to SU(N) 

All of the arguments above are readily extended to theories with an SU(N) gauge group; 
we will only indicate briefly how this goes. The probability distribution Pc{g) i n eq. 
(H) generalizes in the obvious way, with the sum over j replaced by a sum over SU(N) 
representations. Writing 

XF (g) = Ae* (85) 
where F denotes the fundamental representation, A > is real, and G [0, 2ir), let 



n(g) 
2(g) 



int 



Ncf> 



2tt 

exp [27rm (#) / N] 



where int(a;) denotes the integer part of the real number x. This definition assigns a center 
element to every group element, with the property that Z{zg) = zZ(g) for z G Z^. Then 



P c {z)=<<f>[z,U{C)}> where §[z, U(C)] 



1 iiz = Z[U{C)} 
otherwise 



Vt) 
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gives the probability that Z[U(C)] = z. 

Arguments entirely analogous to those in sections 2 and 3 show that the holonomy prob- 
ability Pc{g) approaches the random distribution with only a perimeter-law falloff, while 
Pc( z ) approaches the random distribution with an area-law falloff. At strong-couplings, 
the center elements Z[U(Ci)} fluctuate independently in sub- areas of a large loop, while 
class functions F[U(Ci)] in general do not. From this we conclude that there is magnetic 
disorder among the center elements, but not in the coset. Once again, it should be noted 
that the correlation among SU(N)/Z^ coset elements relies on non-planar contributions, 
which are dominant for large loops. If we would take the large-N limit before the large- 
loop limit, then there is Casimir scaling as in D=2 dimensions, and disorder throughout 
the group manifold. 

Introducing an adjoint Higgs field in D = 3 dimensions, and fixing to unitary gauge, 
we can define the loop observables (in continuum notation), invariant under the remnant 

SUD g r0U p 



where H\ denotes the i-th generator of the Cartan subalgebra in representation r of the 
SU(N) group, and (kk) is just an element of the dim(r) xdim(r) diagonal matrix A(C) 
in this represention. The center elements S[A(C)] are extracted, as above, from the 
phase of xf[A(C)]. 

In the monopole Coulomb gas picture, disregarding the effects of the charged bosons, the 
string tension of < rk\A(C)\rk > depends on both representation r and choice of diagonal 
matrix element (kk). Allowing, however, for the effects of the charged W-bosons, these 
string tensions can depend only on the N-ality of representation r, and are independent of 
(kk). Asymptotically there is disorder in the ZL4(C)] elements, but not in the full ^/(l)^ 1 
group manifold. The distribution of abelian magnetic flux, in the Z N disorder regime, is 
not that of a monopole Coulomb gas. 

5 Discussion 

We have stressed in this article the fact that, while massive virtual particles are often 
irrelevant to vacuum fluctuations in the far-infrared, this is not the case for massive charged 
particles in a confining theory. The screening of external charged sources by quanta of the 
matter field is, of course, a rather trivial point, and allows us to conclude that certain loop 
operators have a perimeter falloff. What may be slightly less obvious is the fact that such 
perimeter falloffs have implications for the probability distribution of large-scale vacuum 
fluctuations also in the absence of external charges. This point is best appreciated in, e.g., 
the 3D Georgi-Glashow model, by imagining an integration, in unitary gauge, over the 
W-bosons and Higgs field, to leave an effective action S e ff[A] involving only the photon 
field. There are no longer any explicit, electrically charged fields left in the action to screen 
multiply-charged abelian loops. Instead, the effect of the virtual W-particles has gone into 




8. 
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altering the Boltzman distribution for vacuum fluctuations of the abelian field, such that 
those abelian configurations which would lead to an area law for even-charged loops in the 
Z2 regime have been suppressed. Thus the effective action S e ff[A] is not only quantitatively 
but also qualitatively different, at large scales, from the QED3 action with a lattice cutoff, 
and a monopole Coulomb gas picture is not adequate to describe the confining vacuum in 
the Z 2 regime. 

The picture we are led to, for the onset of Z 2 disorder in the 3D Georgi-Glashow 
model, is indicated schematically in Fig. ^|. For fixed j3c an d sufficiently large /3r, no phase 
transition is encountered as (3h varies from (3h = (YM 3 ) to fls = 00 (QED 3 ) fl8| . The 
curved solid line, however, represents the breaking of the adjoint string, and the loss of 
"Casimir scaling," while the solid line tailing off in a dashed line represents the breaking 
of the flux tube between double-charged abelian sources. The dashed line indicates the 
complete breakdown of the Coulomb gas picture in unitary gauge, as (3h — > 0. All abelian 
Wilson loops vanish in this gauge in the (3h = limit, although it may still be possible to 
define the n = 2 abelian charge screening distance by extrapolation from non-zero (3h. It 
would be interesting to know where the dashed line terminates. 



OO 



Compact QED 



PQ 




adjoint breaks 



Casimir Scaling 



Z2 Disorder 

Pure Yang-Mills 



OO 



Distance 

Figure 2: Confining disorder, extracted from U(l) and SU(2) holonomies, in the 3D Georgi- 
Glashow model. 



The absence of an adjoint string tension at any length scale, for sufficiently large f3n 
at fixed (PcPr), has been seen in numerical simulations of GG 3 [14], and is easily 
understood. A j = 1 representation quark consists of two components (m = ±1) which are 
double-charged under the U(l) subgroup, and one component (m = 0) which is neutral. 
When the confining field is essentially abelian, the neutral component is dominant, and the 
adjoint loop has no area-law falloff at any length scale. In fact, this gives us an interesting 
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criterion for U(l) disorder in an SU(2) gauge theory, regardless of whether the adjoint 
scalar is elementary or composite. It is required that in the U(l) regime 

1. Even-charged loops have area- law falloff. Otherwise, as we have seen, the loops are 
probing Z 2 disorder. 

2. Adjoint loops have perimeter-law falloff. If not, then abelian neutral components are 
also subject to a confining force, and there is disorder over the entire group manifold; 
not only in a U(l) or Z 2 subgroup. 

In D=4 dimensions, the maximal abelian gauge has been studied extensively in pure Yang- 
Mills theory. This gauge defines a composite adjoint Higgs field, U(l) holonomies, and 
monopole currents. The hope is that confining disorder is U(l) disorder which can be 
attributed, as in GG3, to monopoles. Numerically, however, although double-charged loops 
defined in this formulation have an area-law falloff (cf. ref. ||24j| ) at the length scales probed 
by Monte Carlo simulations, this is also true of the adjoint loops at the same distance 
scales; the second criterion above is not satisfied. 

Returning to the 3D Georgi-Glashow model at large /3#, it is interesting to consider how 
the confining abelian fields are arranged at large scales, where there is Z 2 disorder. It is 
useful to think in terms of the effective abelian theory in fl3~7|), obtained from GG% in unitary 
gauge by integrating out the W and Higg fields. Sef/LA] and S[U, 0] are of course equivalent, 
in unitary gauge, so far as the vacuum distribution of the A-field is concerned. S e ff[A], like 
the GG3 model from which it is obtained, will have instanton solutions corresponding to 
monopoles. However, since the monopole Coulomb gas picture breaks down at the onset of 
Z 2 disorder, it must be that the interactions among monopoles are not really Coulombic at 
long distances, and neither is the field distribution of the corresponding magnetic flux. This 
raises the interesting (although at this stage speculative) question of how the abelian flux 
from monopoles is actually organized, on distance scales characteristic of the Z 2 regime. 

As it is only the sign of cos 0(C) which is disordered in the Z 2 regime, while the effective 
action [A] only involves an abelian gauge field, there is a strong implication that the 
magnetic flux due to monopoles is collimated, at sufficiently large distance scales, in units 
of = Collimated flux of these units, with a stochastic distribution of such "fluxons" 
across the minimal area of a large loop, affects only the sign of odd charged loops, leading 
to an identical string tension for all odd-charged loops, and yielding zero string tension for 
all even-charged loops. This is the proper result in the Z 2 disorder regime. If magnetic flux 
is, in fact, collimated in this way, then a Z 2 vortex picture in this regime is quite natural. 

A particular example of ±27r monopole flux organized into Z 2 vortex configurations of 
±7r flux is shown in Fig. |3[ This is by no means the only possibility. In fact, at large 
the scale L at which monopole flux should be collimated in units of ±7r is actually much 
greater than the average monopole separation R, as seen by comparing eqs. ( [72] ) and (|T4|). 
The illustration in Fig. || might be relevant at lower /3h, approaching the pure Yang-Mills 
limit, when L/R is 0(1). As (3h — » 00, the width of the vortex regions would diverge to 

4 For a discussion of such configurations, in the context of the maximal abelian gauge in YM4, see 
ref. (§]. 
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Figure 3: An example of monopole-antimonopole magnetic flux organized into Z2 vortices. 



infinity, and the monopole Coulomb gas picture is valid at all distances. As (5h is reduced, 
the vortex width decreases. A very attractive possibility is that abelian vortices in GG3 
smoothly transform into center vortices of the pure Yang-Mills theory as (3h —* 0, making 
contact with the ideas of refs. [PHI, and the numerics of refs. [fL0|-T2] . 
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